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ABSTRACT 

A  new  approach  for  the  clear  an  velocity  estimation  in  weather 
radar  is  presented.  A  combination  of  nmiparanietiic  with 
parametric  speetial  analyses  allows  us  to  identity  and  extract 
multiple  processes  caused  by  different  scalterer  types  within  a 
single  radar  resolution  volume.  An  example  of  clear  air  observed 
using  an  S-band  dual  polarisation  radar  is  presented.  Heretofore, 
migrating  birds  and  wind-blown  insects  that  are  mixed  within  each 
resolution  volume  caused  such  data  to  be  unusable  lor 
meteorological  interpretation.  In  this  paper,  we  construct  power 
spectral  densities  ot  polanmetric  variables.  We  me  the  polarimetric 
spectral  densities  to  differentiate  live  scattcrci  types  within  the 
observed  radar  resolution  volume.  We  demonstrate  how  our 
combination  of  non -parametric  and  parametric  spectral  analysis 
can  be  used  to  retrieve  the  true  wind  velocity  in  situations  with 
severe  contamination  by  biological  scalterers. 

Index  Term* — -  Doppler  radar.  Spectral  analysis.  Velocity 
measurement. 

I.  INTRODUCTION 


carried  by  the  wind.  We  demonstrate  these  spectral  techniques  on  a 
clear-air  ease  which  occurred  during  fall  of  21X14  at  the  time  of 
nocturnal  bird  migration. 

2.  RADAR  DATA  AND  WEATHER  CONDITIONS 

Time-series  data  were  collected  with  tlx:  NOAA'NSSL  research  S- 
band  radat  (KOUN)  on  September  7.  2004  at  1 1  pm  local  time  (04 
IT).  The  radar  was  in  dual  potaiizalhw  mode,  simultaneously 
transmitting  and  receiving  waves  of  horizontal  (II)  and  vertical  (V) 
polarizations  at  I  HO"  azimuth  and  0.5''  elevation.  The  pulse 
lepelition  lime  was  7X0  |is  and  the  number  of  samples  available  for 
spectral  analysis  was  128.  Consequently  the  unambiguous  range  Ra 
and  velocity  i'„  were  117  km  and  35  m  s  .  F'air  weather  and  a 
North-North-East  wind  at  5  to  10  m  s  '  were  recorded  on  this  day. 
Howev  er,  the  velocities  registered  by  the  radar  ate  reaching  30  m 
s  1  I1M2).(9).  Radar  meteorologists  attribute  this  inconsistency  to 
"contamination  by  biological  scalterers"  aiul  discard  such  data  as 
worthless  for  meteorological  interpretation.  We  demonstrate  a 
technique  that  allows  us  to  retrieve  both  the  wind  and  bird  speeds 
in  this  and  simitai  sitmliuns. 


The  national  network  of  Doppler  radars  (WSR-88D)  estimates  the 
vertical  profiles  of  wind.  These  profiles  are  obtained  whenever 
back-scattering  is  sufficient  to  produce  delectable  signals.  In  the 
alncncc  of  clouds,  such  backscaltcnng  is  prov  ided  by  insects, 
irregulanties  in  icfradivc  index,  and  smoke.  Small  insects  could 
potentially  bias  wind  velocity  estimates,  though  this  bias  is 
generally  considered  to  be  insignificant.  limlv.  however,  introduce 
large  biases  to  the  Doppler  wind  measurements  especially  during 
migration  (3|,(6).(9).(H)).|I2|.  Radar  engineers  continually  try  to 
improve  the  resolution  of  radar  products  |4).(S),(II|.(I2|.|I3). 
When  there  arc  two  ot  more  processes  ongoing  in  a  resolution 
volume,  these  might  not  he  recognized  in  the  finest  resolution. 
Migrating  bird*  bias  weather  velocity  and  reflectivity  estimates 
brutally,  often  requiring  users  to  censor  or  discard  the  radar  data. 
We  analyze  Doppler  spectra  and  compose  fields  of  spectral 
densities  of  polanmetiic  variables  which  expose  the  potential  for 
obtaining  measurements  of  both  insects  and  hitds.  Consequently 
the  insects'  measurements  lead  to  ohtauimg  winds  aloft,  and  the 
bird  measurements  provide  inhumation  on  the  bird  activities.  We 
use  autoregressive  models  to  demonstrate  our  potential  for 
identifying  tlx*  many  types  of  sealterets  with  diflerciu  velocities 
within  the  resolution  volume.  This  opens  possibilities  to  not  only 
separate  birds  from  insects,  but  also  to  separate  two  types  of 
insects,  e.g.  the  insects  biasing  the  wind  estimates  and  the  insects 


3.  POLARIMETRIC  SPECTRAL  DENSITIES 

The  power  spectral  density  (PSD)  is  estimated  in  two  ways.  First, 
the  PSD  is  estimated  from  tlx-  DFT  of  the  time  series  data 
weighted  with  the  von  Hami  (Hanning)  window,  tlx-  ground  clutter 
is  removed  with  a  "notch  and  interpolate"  filter  at  frequencies 
corresponding  to  velocities  between  -I  in  s  1  and  I  ms  1  [5).  The 
spectral  noise  level  is  estimated  using  rank  order  statistics  on  tlx- 
speCral  coelliex-nts  at  every  tange  location  and  then  removed. 
Second  the  echo  signal  is  found  from  the  inverse  DFT  of  tlx- 
clutter  filtered  PSD  and  used  lor  parametric  PSD  estimation.  The 
paiamctnc  PSDs  are  approximated  with  the  multi-signal  spectral 
decomposition  technique  using  the  autoregressive  models  ARfp) 
with  p  -  2.  3.  4  (5).(7).  In  a  dual  polarization  system  both  tlx-  H 
and  V  polarization  spectral  densities  are  available;  from  these  it  is 
possible  to  estimate  the  differential  properties  at  localized  resolved 
Doppler  shills  | l|,(2J.  For  .W  transmitted  pulses  there  are  M 
spectral  cocfticx-nts  in  the  spectra  from  a  resolution  volume. 
Differential  reflectivity  Znr.  between  II  and  V  channels  is  defined 
as  a  diflcrence  (ui  logarithmic  scale)  of  the  mean  powers  in  these 
channels.  To  obtain  spectral  densities  of  dilTercntial  reflectivity, 
the  values  are  computed  for  every  H-V  pan  of  spectral 
coefficients  of  the  power  spectral  densities  |2).|X].|  10|.  The 
spectral  densities  of  the  complex  co-polar  correlation  coefficient 
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{magnitude  and  ptusc  jnr)  are  estimated  from  a  running  3- 
point  average  on  contiguous  complex  spectral  coefficients  of  the  II 
and  V  channel  |2J.  The  backscattcr  diffcrcntul  phisc  Sit  found  by 
subtracting  ihc  system  phase  from  The  system  phase  of  the 
radar  digital  receiver  is  found  from  the  ground  clutter  returns  (15]. 

4.  MONPARAMETRIC  spectral  densities  in  range 

f  ields  of  Doppler  spectral  densities  arc  presented  as  a  function  of 
ranee.  The  Nooparamciric  PSD  in  Ihc  II  channel  (S.)  and 
polarimctric  spectral  densities  (ZM  .  and  .*»  alcmg  the  ISO1 
radial  are  shown  in  Fig.  la.  b.  c.  d.  Only  spectral  coefficients  with 

5,  '  5dll  above  the  noise  are  displayed.  One  column  m  the  image 
at  any  range  represents  the  spectrum  at  that  range.  In  the  presence 
of  scaltcrers  the  image  is  expected  to  display  a  somewhat 
continuous  distinct  hind  corresponding  to  the  radial  component  of 
the  scattercrs'  velocities.  Thus,  in  an  ideal  case  this  band  represents 
the  dependence  of  radial  velocities  on  range  (height  I  that  should  be 
consistent  with  the  soundmg  wind  profile.  Other  curves  and  blobs, 
which  deviate  from  the  band  tracing  the  mean  wind,  arc 
contaminants  (l|.  The  2 -dimensional  histograms  of  S,.  /fix.  <x- 
and  S  (Fig.  Ic.  f.  g.  h)  constructed  for  ranges  50  to  SO  km  reveal 
the  distribution  of  each  parameter  in  Ihc  velocity  space.  The 
density  fields  disclose  two  bands  located  at  about  10  ami  20  m  s1 
respectively. 

The  band  at  10  m  s'  is  caused  by  insects.  The  insect  band 
appears  as  a  dark  wide  hand  in  the  spectral  fields  and  shows  up  as 
a  confined  blob  m  the  polarimctric  histograms.  The  insect  bard 
exhibits  smaller  S,  (up  to  25  dB);  higher  value*  of  ZM  (3  <ZM  *10 
dB)  and  values  above  0.8  with  most  occurrence*  close  to  I .  ard 
a  narrow  distribution  of  phases  |40:  *  S<  80°  |. 

The  hand  at  20  m  s'  rs  caused  hy  birds  The  bird  hand  can  he 
characterized  by  larger  powers  (S,  *  41MB).  smaller  10  <ZM  < 
8  dB).  wide  distribution  of  />.  (0.5  <  <  Ik  and  all  possible 

pluses  (  ISff  <  S<  1 80°)  PolarimelrK  propertses.  shapes  of  the 
hands  and  Ihc  values  of  the  associated  velocities  reveal  that  insects 
and  birds  shared  the  airspace  throughout  most  of  th:  boundary 
layer  on  that  night  (2|.  The  polarimctric  signatures  of  insects  ard 
birds  arc  substantially  different  (12],  therefore  we  can  discriminate 
the  two  types  of  scaltcrers  if  they  are  separated  in  vclceily. 

5.  PARAMETRIC  SPECTRVL  DENSITIES 

In  weather  radar  two  spectral  charactcrrstics  arc  of  main 
importance:  th:  power  for  the  reflectivity  and  the  peak  frequency 
for  the  vclceily  computations.  Paramclnc  spectral  techniques  such 
as  Burg,  eigenvector  decomposition,  MUSIC  I  multi  signal 
classification),  etc.,  distort  power  content  but  can  provide  higher 
resolution  frequency  cslitmtes.  The  MUSIC  algorithm  computes 
the  frequency  spectrum  of  a  signal  by  decomposing  it  in  complex 
exponentials  superimposed  cm  noise  (7|. 

The  mean  velocitars  of  the  hind*  with  the  diverse  scattercrs  are 
visible  in  the  plots  of  spectral  fields  (Fig.la.  let  and  can  he 
estimated  using  either  spectral  coefficients  (I]  or  spectral  VAD 
(velocity  a/imuth  display  I  analysis  (2).  We  use  a  different 
approach  here  and  estimate  these  velocities  using  parametric 
spectral  modeling.  Two  processes  are  obvious  in  the  clutter  filtered 
data  set:  th:  wind  traced  by  insects  and  Ihc  migrating  bind  flow’. 
The  scalKrcrx  of  the  same  type  tend  to  have  similar  velocitKs. 
Since  we  can  guess  the  number  of  processes,  we  can  model  the 


signal.  Thus,  two  sinusoids  capture  mean  speeds  for  birds  and 
msec!*,  and  three  sinutoeds  can  reveal  additional  passive  insects.  A 
spectrum  at  .<0  km  range  (fig.2a)  consists  of  three  portions  from 
biological  scattercrs  (Fig.2c.  d.  c)  bads  and  two  types  of  insects. 
Thus  it  can  be  modeled  with  3  sinusoids  as  shsiwn  in  l  ie— h. 

SlA) 


A.,  ""i 


A.  "HI 


Fig.  I  Spectral  drnsities  of  a)  Sj.  bl  Z,„.  c)  ,\ ..  d)  A  and  their 
histograms  c)  Si.  0  <f.v.  g)  /*,.  and  h)  S .  Only  spectral  cceflicKnts 
with  S.  >  5dB  above  the  noise  are  displayed.  Bads  signature*  show 
-20  m  s'1  velocity,  smaller  powers  <S,  *  25  dB|:  smaller  ZM  (0  < 
Zoe  <  8)  dB.  kiwer  /*,  (0.5  *  p>,  <  I).  and  all  possible  phases  ( 
ISO1  <  S<  ISO”).  Insect  signatures  show  -10  m  s'1  velocity,  larger 
powers  <&  *  41)  dB):  larger  Zb,  (3  <  Zt*<  10  dB  l.  higher  <0.8 
*  1).  and  a  narrow  span  of  phases  (40:  *  S<  80”). 


l  it!.  2.  Spectral  densities  from  various  processes  in  the  resolution 
volume  a|  a  composite  spectrum  at  30  km  tFig.l).  b)  spectral  lines 
representing  complex  cxporcntials  in  noise  corresponding  to  the 
mean  values  of  composites:  e)  portion  due  to  migralotg  birds:  d) 
portion  due  to  active  insects:  cl  portion  die  to  wind  tracers. 
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The  half  of  Ihc  PSD  from  Fig  la  with  positive  velocities  it 
repealed  in  Fig. 3a  fee  visual  assessment  the  PSDs  estimated  using 
autoregressive  models  of  order  2.  3.  and  4  are  presented  in  Fig  3b. 
c.  and  d  respectively.  The  ARI2I  model  shows  two  tfc finite  bands 
front  which  the  velocities  of  the  two  sealtercr  types  can  be  easily 
estimated.  The  AR(3|  and  AR|4|  mixfcls  also  show  2  bands  but  the 
velocities  corresponding  to  these  band  ate  less  prominent  Fig. 4 
shows  zoomed  in  pontons  of  Ihc  non-parametric.  AR)2)  and  AR(3) 
PSDs  The  AR|3)  model  appears  to  split  the  insect  band  into  two 
bunds  at  about  5  m  s'  and  10  lit  s'1.  The  atmospheric  sounding 
(rawinsoode)  at  0  UT  detected  7  m  s'1  wind  at  the  surface  |2|. 
However,  the  velocity  estimated  using  l be  spectral  VAD  technique 
at  4  UT  and  541  km  range  was  12ms'  in  the  direction  of  the  wind 
(220°)  and  9  m  s  1  m  the  ISO0  radial  [2).  Hither  the  wind  increased 
from  0  UT  to  4  UT  or  Ihc  insects  were  active  flyer*.  The  AR<2> 
estimate  has  a  peak  at  7  m  s’  and  341  km  tango.  ARI3)  medels 
approximate  two  peaks  (at  5  and  10  m  s'1)  at  341  km  range.  These 
appear  to  be  the  trace  of  passive  insects  carried  hy  the  w  ind  and  the 
trace  of  the  mean  velocity  of  active  insects. 

The  Ziu  between  the  mixfclcd  channels  (Fig.  5)  exhibits  lie 
polarimctric  ditlercnces  in  hands  with  higher  values  for  Ihc  insect 
bund.  Although  it  loc&s  correct  in  lie  presented  example,  we  do 
not  think  su;h  a  Zm  is  trustworthy.  If  the  peaks  in  the  modeled 
PSDs  of  II  and  V  channels  do  not  exactly  match  then  the  7.m 
values  may  vary  dramatically.  We  suspect  that  the  AR(2l  model 


capture*  live  wind  and  bird*  nicotic*  accurately,  a*  determined  by 
independent  observation*.  It  teem*  that  the  AR(J)  model  ccwrcctly 
obtain*  the  velocity  of  the  wold-carried  insects.  These  are  the 
insects  that  hia*  the  wind  velocity  by  2  to  3  m  s'.  Ihe  AR(3| 
model  also  appear*  to  accurately  estimate  the  herds  velocities. 
Evidently.  the  AR(4)  model  has  too  high  of  an  order,  and  thin 
captures  spurious  details  instead  of  the  mean  speed  of  the  different 
scatterer*.  We  submit  that  the  combination  of  the  AR(2)  and  AR(3| 
model*  arc  sufficient  for  wind  velocity  estimation  in  the  frequency 
domain  m  ca*e*  with  severe  contamination  by  migrating  birds  and 
similar  situations. 

6.  PROCESSING  TECHNIQUE 

The  reason  that  sve  use  a  combination  of  non-parametric  and 
parametric  ipcctral  analyses  i*  bused  on  the  following.  I:ir*t.  the 
non-parametric  technique  is  well  known  and  accepted  for  the 
weather  radar  application*  |5).  Second,  the  parametric  technique 
we  use  require*  a  prior  knowledge  of  the  physical  process  [7) 
which  the  non-puranvetric  method  can  provide.  Third, 
meteorological  parameter*  *u:h  a*  spread  of  Doppler  velccitK*. 
and  reflectivity  arc  difficult  to  obtain  with  parametric  methods. 
And  last,  at  low  elevation*  the  ground  clutter  is  sporadic  and 
therefore  the  model  order  for  matching  the  physical  process  it  not 
known  unless  there  is  a  priori  or  spectral  information. 

The  presented  technique  is  sumnun/ed  below: 

1.  Intimate  non  parametric  .spectrum  (DFT)  fee  each  restitution 
volume. 

2.  Apply  the  Smirnov- Kolmogorov  test  on  the  spectrum 
coclTicicnls  to  estimate  the  noise  power.  Then  remove  the  none. 

3.  Check  power  of  spectral  coefficients  at  and  close  to  zero 
velocity;  notch  them  if  gmmd  clutter  is  detected:  interpolate  across 
Ihe  notch  using  adjacent  non-clutter  cocA~icic«s. 

4.  Compute  pobnmctnc  spectral  densities  from  the  spectra  of  the 
horizontal  and  vertical  channels  [2|. 

5.  From  pokuimclnc  spectral  densities  estimate  the  number  of  the 
scatterer  types 

6.  Perform  the  1DFT  on  the  grmmd-cluttcr  corrected  signal. 

7.  Estimate  spectrum  using  the  parametric  model  of  the  order 
estimated  in  5. 

8.  Extract  the  wind  velocity  from  the  hand  of  passive  scatterer*  in 
th:  spectra  using  a  special  technique,  such  as  spectral  VAD  (2). 

7.  CONCLUSION 

Polarimelric  spectral  analysis  of  the  w  eather  radar  returns  in  clear 
air  and  at  night  reveals  the  simultaneous  presence  of  migrating 
birds  and  msecls  mixed  within  radar  resolution  volumes.  The 
spectral  signatures  of  polarimelric  variables  can  identify  the 
scatterer  types  while  the  power  spectral  densities  estimated  with 
autoregressive  models  provide  the  velocity  of  the  different 
scatterer  types.  The  wind  velocity  can  be  estimated  from  the  wind 
hand  corresponding  to  the  insects. 

Wc  lave  demonstrated  that  an  autoregressive  model  can 
pinpoint  several  velocities  of  ccmlrlbuting  scatterer*  within  the 
sain:  resolution  volume.  We  believe  this  rt  the  first  attempt  to 
identify  velocities  of  three  (not  cmc)  types  of  contributing  scatterer* 
for  a  single  resolution  volume.  We  have  shown  that  combination  of 
polarimelric  non- parametric  with  parametric  spectral  analyses  can 
identify  birds,  insects  and  wind. 
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